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Abstract 

The paper presents an evolutionary game-theoretic approach to open 
access publishing as an asymmetric game between scientists and pub- 
lishers. We show how the ordinary differential equations of the model 
presented can be written as a system of Hamiltonian partial differential 
equations. The understanding of the setting as a Hamiltonian system im- 
plies some properties reflecting the qualitative behavior of the system. 

1 Introduction 

The topic of open access publishing has been extensively and controversially dis- 
cussed. For general information see, for instance, [TTl [121 US] . There are many 
models of open access named after different colors such as golden or green roads 
to open access [21 [10], but we do not want to go into more details, here. Con- 
cerning scientific publishing particularly in mathematics, wc refer to [2 [7] and 
[S] as well as to the references therein. 

In [4] the authors illustrate a game-theoretic approach to open access publishing 
in order to understand different publication patterns within different scientific 
disciplines. The underlying observation is that there are communities where 
open access publishing is widely adopted, whereas other scientific communities 
are far away from practicing any open access publishing. First, different classical 
game settings are discussed namely a zero sum game, a game similar to the 
Prisoners' Dilemma (up to sign), and a stag hunt game version, all describing a 
Nash equilibrium dilemma of the non-open access communities. Second, these 
classical settings are transferred into their quantum game extensions which allow 
to tackle the dilemma mentioned that cannot be solved within the classical 
setting. 

The classical approach in [4] is formulated as a two-player game where the 
players are authors (scientists) in symmetrical situations. The two authors have 
the same set {si, S2} of strategies at their disposal. Consequently, the game looks 
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exactly the same to both of its players. For the open access game each player 
has to decide between the option to pubhsh open access as the first strategy si 
or, as the second strategy S2, to conventionally publish in traditional journals 
where articles go through peer reviewing. Depending on the strategy chosen by 
the co-player, the players aim to maximize their success of the game. Due to 
the symmetric situation both players do not only have the same set of strategies 
but also have the same payoff matrices. 

In contrast to the approach of a symmetric two-scientists game, the present 
paper attempts to explain mathematically the open access play as a conflict 
of interest between scientists and publishers, as it is emphasized by the rapidly 
accumulating literature on open access and electronic publishing. Assuming that 
the open access problem is not a conflict between those scientists who publish 
open access and those who do not, one takes the view that it is more adequate 
to understand the problem of open access as an asymmetric conflict between 
scientists and publishers. 

The crucial point of this view is that it actually suggests a bimatrix game 
describing this asymmetric setting. 

We do not explain why both players - scientists and publishers - behave as they 
do in the game situation described, but we will give a mathematical description 
of their interaction in the game. Moreover, we do not argue for preferences in 
choosing one game strategy over another. 



We consider scientists and publishers as players in different positions them hav- 
ing different strategy sets as well as different payoff matrices. Moreover, we take 
into consideration whole populations of players parts of which choose either one 
or the other strategy. 

The set of pure strategies for the scientists is {si, S2}, where si is publishing open 
access and S2 stands for conventional publishing. Furthermore, let < xi < 1 
and < a;2 < 1 denote the relative parts of the scientists' population playing 
strategy si and S2, respectively. If, for example, half of all scientists publishes 

open access and the other does not, one has xi ~ X2 ~ —■ Or. if the relation 
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is 20% to 80%, then xi = — and 2:2 = -. Since xi and X2 describe relative 
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frequencies, we see that xi + X2 — I always holds true. Hence, the set of all 
possible mixed strategies for the scientists' population is 



The opponent population of publishers consists of two types playing pure strate- 
gies {pi,p2}. Here strategy pi means the publisher accepts or realizes open access 
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publishing, whereas p2 represents the strategy of definite dechning any way of 
open access. For the publishers' population let yi be the frequency of strategy pi, 
where i = 1,2. Thus, all possible mixed strategies for the publishers' population 
are given by the set 

^ = |y = (^y^ e M2 : < yuy2 < l and yi + = l| • 

In order to discuss the conflict of open access publishing, we proceed with estab- 
lishing the payoff matrices A for the scientists and B for the publishers. Each 
payoff depends on the strategies chosen by the participating players. A player 
in the scientists' population using strategy Si against a player from the publish- 
ers' population using strategy pj obtains the payoff a^ , whereas the opponent 
obtains bji, with i,j = 1,2. 

As reasoned in [4], it is convenient to assume that the scientist tries to maxi- 
mize his scientific reputation. Let R > denote the reputation payoff that the 
scientist is awarded. This payoff will be slightly reduced by < r < i? in case of 
publishing open access. Hence, we still have i? — r > for the total reputation. 
In addition, we take into consideration the impact of an article by the scientist. 
The impact gives some payoff / > both for the scientist and for the publisher. 
Further, this impact obviously will reduce by some < t < /, if the concerning 
journal is not sufficiently available. Moreover, we assume that l < r saying that 
the reduced impact in case of non-open access publishing is smaller than the 
loss of reputation in case of open access publishing. Open access publishing will 
cause expenses i > that are shared equally by both players if both play their 
open access strategy. Furthermore, G > expresses some moderate journal price 
scientists have to pay for library subscriptions as well as the compensation for 
expenses of the publisher. Finally, let us introduce P > for exorbitant profit 
representing the prices of the most expensive journals coupled with continuing 
non-open access publishing and taking the tremendously increasing fees for ac- 
cess for granted. For our considerations, we make the realistic assumption that 
P is so large that the inequality 

G + P - L> r - i (2.1) 

holds true, which is equivalent to 

G + P> L + r-L. 

We already have r — i > and the open access expenses L actually cannot 
exceed the upper bound given by this inequality. 

If an open access publishing scientist meets a publisher accepting open access, 
the scientist will get his reduced reputation, his research impact, and will pay 
for maintaining some open access publishing system as well as for the journal 
subscription. The publisher takes the price for the journal subscription and some 
amount given by the impact. Apart from this, the publisher contributes to some 
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extent to open access. Hence, the payoff is (i? — ?■)+/— — — G for the scientist 
and G + I — — for the pubhsher. 

If the scientist's strategy is open access whereas the pubhsher plays the non- 
open access strategy, the scientist's payoff wiU be (i? — r) + / — L. The scientists 
have to run the open access completely by themselves and do not have any 
expenses for journal subscriptions. For the publisher nothing happens and hence 
the payoff is 0. 

A non-open access publishing scientist gets his full reputation R but only re- 
duced research impact / — t. Additionally, the journal has to be paid, which 
reduces the payoff by G. If the publisher's counter-strategy is open access, 
the payoff for the scientist will be R+ {I — l) — G. In this case the payoff 
for the publisher will be G + {I — i) — L paying the total sum for open access 
acceptance. Moreover, the payoff is R+{I — l) — G — P for the scientist and 
G -I- (/ — i) + P for the publisher in case of the non-open access strategy of the 
publisher. 

Altogether, we have the complete income-and-loss statement as given in Table[TJ 



Table I: Complete payoffs 



strategies 


payoff 
for the scientist 


payoff 
for the publisher 


Sl < — > pi 


{R-r)+I-^-G 




Sl < > P2 


{R-r) + I -L 





S2 < > Pi 


R+{I-l)-G 


G+{I-l)-L 


32 < ^ P2 


R+{I - l) - G - P 


G + {I-i) + P 



Hence, the corresponding payoff matrices are 



A 



, {R-r)+I-^-G (R-r) + I-L 

ail ai2 \ _ \ ' 2 ^ ' 



0,21 CL22 



R+{I -l)-G R+{I-l)-G-P 
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and 

B 



On bi2 \ _ \ 2 

^21 ^22 






Now, both payoff matrices determine tfie payoffs for the entire populations. For 
a general description of how this happens we warmly recommend the book [S] 
by Hofbaucr and Sigmund. The corresponding strategies are spread within both 
populations in accordance with the frequencies of the strategies. Hence, if the 
scientists' population is in state x Cz S and the publishers' population is in state 
y € P, then the payoff for the entire population of scientists will be 

{x,Ay) = x^Ay 

and that for the publishers' population will be 

{y,Bx) = y^Bx = x'^B^y . 



As already mentioned, players will choose their game strategies with the inten- 
tion of maximizing their average payoff. This would be easy if one knew the 
strategy that the opponent player is going to choose. In case both players of the 
game choose simultaneously strategies of best reply to the choices of the others, 
a pair of strategies occurs where both players are incited to keep going into the 
same direction. Such a pair of strategy choices is known as Nash equilibrium. 
This is technically expressed by the following definition. 

Definition 2.1 A Nash equilibrium is a pair {x,y)ESxP such that 

Ay < x^Ay for all x e S 

as well as 

y^Bx < y^Bx for all y e P 

is fulfilled. 

The following considerations aim to give a Nash equilibrium for our open access 
game setting. 

Lemma 2.2 Each of the differences ai2 ~ 022 • 121 ~ an, ^22 — &21, andbu—bri 
is positive. 

Proof. Indeed, 

ai2 - a22 = - L + L + G + P > 

by (El]) and 

L 

021 - an = -t + r + — > 
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by r > t. Furthermore, 

b22 - = G + I - L + P 
which is positive since I > t. FinaUy, 

611 - bi2 = ^ + L 

which is obviously positive. □ 

Lemma 2.3 Let {x,y)<ESxP be given by 

- f xo \ , _ f ijo 

x = [ and y = k 



whe- 



re 



622 — bi2 ai2 — a22 

^0 - 7 r — rr J— ('^d yo 



^22 ~ bi2 + ^11 ~ ^21 0-12 ~ 022 + 121 ~ All 

Then 

x^Ay = x^Ay for all x e S 

and 

y^Bx = y^ Bx for all y £ P 

hold true. 



Proof. First, by Lemma [272] both denominators of the fractions defining xq and 
yo are positive. Then we have 



Ay 

as well as 
Bx 



an ai2 \ I yo 
0.21 0,22 J \1 — yo 
yo(aii - ai2) + ai2A det A 



2/o(a2i - a22) + a22/ ai2 - a22 + 021 - an yl 



bii 612 \ f Xo 
b21 b22 J \l - Xo 
xoibn - bi2) + bi2\ - det B fl 



xaib2i - ^22) + ^22/ 622 - bi2 + bii - &21 yl 



Hence, 



rp , det A , , , , det A 
x^ Ay = (x + (1 - x)) ^ 

ai2 - a22 + a2i - an ai2 - 022 + 021 - an 
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docs not depend on a; = ( I and 



V^Bx = —{y + (1 - y)) = 



dctB , , _ detS 



^12 - ?'22 + ^21 - &11 ' ' br2 - b22 + hi - hi 

is independent of y = ( ^ ) . □ 

Corollary 2.4 The pair {x,y) E S x P given in Lemma \2.3\ is a Nash equilib- 
rium. □ 

Remark 2.5 For this game there is no strict, i.e. given by pure strategies, Nash 
equihbrium. The Nash equiUbrium {x^y) Cz S x P specified in Lemma 12.31 is a 
unique mixed Nash equilibrium and determined by the foUowing equations, c.f. 
[5], Chapter 10.2 (Page 116, above), 

aiij/o + ai2(l - 2;o) = a2iyo + 022(1 - yo) 
611X0 + &i2(l - a^o) = 62ia;o + ^22(1 - 2:0) • 

□ 



3 The Dynamics of the Game 

Due to the payoffs of the game, there is an intrinsic dynamics within the system. 
Actually, starting for instance with the situation that the scientist has chosen 
strategy si whereas the publisher's strategy is pi, the outcome is 

{R-r) + I-^-G 

for the scientist and 

for the publisher. In this situation the scientist could increase his payoff by 
changing his strategy to S2, since r — t > and L > 0. Indeed, the inequality 

R-r + I-^-G<R + I-L-G 

is equivalent to 

L 

which holds true obviously for r — i > and L > 0. After this change in the 
scientist's strategy the payoff is 
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for the scientist. Holding strategy pi the pubhsher's payoff is 

G+iI-L)-L. 

In this situation, the pubhsher sees the chance to change strategy to p2- The 
reason is that the pubhsher's payoff of G + (/ — i) — i in the game S2 <-> pi can 
be increased to G + (/ — /,) + P in a game S2 ^ P2- In fact, L > and P > 
give 

G+{I-i)~L<G^{I-i) + P. 

However, finding oneself in the situation of the game S2 P2 the scientist 
notices that G and P are unnecessary costs, L would not be so huge and the 
loss of reputation will be moderate. Hence, the scientist will prefer to choose 
strategy si. Altogether, changing strategies in order to maximize outcome ends 
up moving in a circle. 

The dynamics for pairs of pure strategies looks like 

(si,Pl) ^ {si,P2) 

i T 
(s2,Pl) -> (S2,P2) 

where the horizontal arrows picture the change of the publishers' strategies and 
the vertical arrows illustrate the deviation in the scientists' strategies. 

This dealing with an "oscillating" system requires a non-static approach. There- 
fore, we consider the frequencies of the strategies within both populations 
as time-dependent quantities which are differentiable functions of i G M, i.e. 
Xi = Xi{t) and j/j = Vji^) for *i J = Ij 2. The first derivatives of these functions 

Xi = — — and Vi = —r- for i, 1 = 1,2 
dt di 'J > 

describe the rate of growth of the respective frequencies. 

Modeling the dynamics and the cyclic behavior of that system, we follow the 
ideas of Hofbauer and Sigmund developed in [5], Chapter 10.3. The increase of 
strategy Sj of the scientists' population is given by the per capita growth rate 

— and equals the difference between its payoff {Ay)i and the average payoff 

Ay. The same holds true for the publishers' population. This leads to the 
system of ordinary differential equations 

Xi = X, {{Ay), - x^Ay) 

for i,j = 1,2. 

Let us look closer at these differential equations. We have 

j\^y ^ f «ii "12 ^ (^^1^ = fanyi+ai2y2\ 
\ 0.21 022 J \y2) \a21y1 + a22y2 J 
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and 



Hence, 



( bn 


b2i \ 






\ bi2 


b22 J 




-{ 



biixi + 621X2 
^122:1 + 6222:2 



iAy)t - x^Ay = yiian - xiau - X2a2i) + 2/2(012 - xiai2 - ^2022) 



and 



{B^x) 



y B a; = a::i(6ij - j/1611 - 1/2^12) + X2(62j - 1/1621 - 2;2&22) 



Since xi + a:2 = 1 and j/i + 2/2 = Ij it is admissible to introduce new variables 
X := Xi and y := j/i (then X2 = 1 — a; and 2/2 = 1 — y) and consider the 
corresponding differential equations for x and y. These are 

X = ±1 

= xi {{Ay)i - x'^Ay) 

= X (2/1 (an - xioii - 0:2021) + 2/2(012 - a::iai2 - a:2a22)) 

= X (2/(011 - xoii - (1 - a::)o2i) + (1 - 2/)(oi2 - 2:012 - (1 - 2:)o22)) 

= x{l - x) (2/(011 - 021) + (1 - 2/)(oi2 - 022)) 

= x{l - x) ((ai2 - 022) - (ai2 - 022 + 021 - 011)2/) 

and similarly 

2/ = 2/1 = 2/(1 - y) (-(^22 - ^21) + (^22 - ^21 + ^11 - ^12)2;) • 

Now, let us introduce new constants by 



a 


■ = ai2 


- 022 


b 


■ = 021 


- an 


c 


:= ^22 


- 621 


d 


:= &11 


- 612 



which are positive by Lemma 
namely 



Hence, only two equations are remaining 



X = x{l — x)(a — {a + b)y) 
y = 2/(1 - 2/)(-c+ (c + d)j:) 



(3.1) 
(3.2) 



This is an adequate form for investigating the qualitative behavior of the system. 

The subsequent section is devoted to a further investigation of the dynamics of 
the system given by Equations (|3.ip and 
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4 Symplectic Reformulation 



We consider the square Q = {(x, y) G : < x,y < 1}. Obviously, Q = SxP 

by 

We use symplectic techniques for investigating what happens in the inte- 
rior intQ of Q. More precisely, in order to prove that the orbits of the Sys- 
tem p.ip and (|3.2[1 are closed and do not leave intQ we describe the whole 
setting as a Hamiltonian system. 

The interior of Q is the open subset 

M = {{x,y) G : 0<x,y <1} 

of M^. The function : A/ — > R given by (p{x, y) = xy(\ — x)[l — y) is positive. 
Hence, 

Lo — —dx A dy 

defines a 2-form a; on M. 

Since we are in dimension 2 the following is clear. 

Lemma 4.1 The pair {M,uj) is a symplectic manifold. □ 

We are going to derive that Equations (|3.ip and p.2p are the equations of motion 
of the Hamiltonian system (M, where the Hamiltonian H : M ^ R is 

defined as 

H{x, y) = cln(a-) + (iln(l - x) + aln(?/) + &ln(l - y) for {x, y) G M . 

The notion of Hamiltonian systems is well-known and widely used in mathe- 
matical physics to describe the behavior of systems in classical mechanics. In 
physics the Hamiltonian H describes the energy of the system under consider- 
ation. Hence, it would be interesting to understand, what intrinsic property is 
determined by the "energy" of open access publishing. 

Lemma 4.2 The Hamiltonian vector field of {M,uj,H) is given by 

( \ 

"Pit 

dy 

Xh = 

dH 
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Proof. Generally, Xh is defined by the equation 

lo{Xh, )=dH . 

Clearly, 

dH ^ dH ^ 
ox oy 

If V is any vector field on M with functions Vi and V2 as components, then 

Uj{V, ) = —dx A dy{V, ) = —{vidy — v-)dx} . 
tp if 

Hence, 

) = dH 

if and only if 

dH I dH I 

-TT- = V2 and -T— = —vi , 

ox (p oy (p 

which proves the assertion. □ 

Let us remind that motions of the Hamiltonian System (Af , oj, H) are curves 
7(i) in M satisfying the differential equation 

Altogether, we now can prove the following proposition. 

Proposition 4.3 The solutions of iS. 1]) and i3.2\) are exactly the motions of 
the Hamiltonian system {M,lu,H). 



Proof. We have 



dH \ a 

if— = xy{l - x){l - y) < - 



and 



dy [y l-y 

ax{l — a;)(l — y) — bxy{l — x) 
x{l — x){a — ay — by) 



-(/?—— = -xy{l - x)(l - y) 



dx IX 1 — X 

= -~cy{l ~ x){l - y) + dxy{l - y) 
~ y{l — y){—c + ex + dx) . 

Thus, Equations (|3.ip and (|3.2p are equivalent to 

dH dH 

ip— = X and - ip— = y , 
dy dx 
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which implies the statement by Lemma 14.21 □ 

As a consequence, the solutions in M of the System (|3.1|1 and (|3.2p correspond 
to level set curves of the Hamiltonian i7, i.e. the Hamiltonian H is constant 
along solutions. 

Corollary 4.4 Let {x{t),y{t)) give a solution of Equations i3.1\) and US. 2]) . 
Then 

^H{x{t),y{t))^0. 



dt 



Proof. Indeed, 



±Hixit),yit)) = 

= ^ {cln{x{t)) + dln(l - x{t)) + aln{y{t)) + 61n(l - y{t))} 

X X y T y 

= c a- h a 0- 



X 1 — X y l^y 
= (c(l - x) - dx)(a - (a + b)y) + (a(l - y) - by){-c + (c + d)x) 
= (c — cx — dx){a — ay — by) + {a — ay — by){—c + cx + dx) 
= 

by using (|3.ip as well as p.2p for replacing i and y. □ 

In the remainder of this section we want to discuss the shape of the level set 
curves of H. In particular, we will show that these curves arc closed, which 
implies that the solutions of (|3.ip and (|3.2p are periodic. 

First, we are going to determine the extremal points of H . By 

dH a b 



and 



we obtain 



dy y 1-y 
dH c d 



dx X 1 — X 



a b c 

and 



yo 1 - yo xo I- xo 

as necessary conditions for a critical point (xo,?/o)- This leads to 



a c 
yo = — —f and Xq = — — 

a + b c + a 



which corresponds exactly to our Nash equilibrium described in Lemma l2.3l By 







c 


d 




d'^x 


x2 (1 






d^H 


a 


b 




d'^y 


y' (1 


-y? 






= 0, 




dxdy 


dydx 





and 
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we have 

(c + d)3 

^(.To,2/o) = ^ < and 

d'^H. , (a + 6)3 

at the point (a;o,j/o)- This impUes that the matrix 

/ {c + df ^ 

V ab 
is negative definite showing that H takes in (xq, yo) a maximum as the unique 
extremum. Second, since 



we have 



Um ln(t) = —oo 



hm H{x, y) — — oo 

{x,y)~,q 



for each q in the boundary dQ of Q. Hence, the graph of the Hamihonian H 
looks hke a hill over the square Q. 




Figure 1: Graph of the Hamiltonian H with parameters a — 1, b — 2, c — 2, and d — 3 



The top of the hill represents the unique maximum. Towards the boundary dQ, 
the hill crashes into bottomless depths. 

Since all solutions of our system correspond to level set curves, they are visual- 
ized as contour lines of the hill described. Hence, all orbits surround the point 
(xq, yo) and remain in M. 
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Figure 2: Periodic orbits surrounding the Nash equilibrium (xq, yo) with parameters a = 1, 

6 = 2, c = 2, and d = 3 



5 Further Discussion 

Obviously, the model presented here is a basic approach to discuss the phe- 
nomenon of multiple arrangements of realizing open access publishing in differ- 
ent scientific communities. However, in contrast to the approach presented in [1] 
it treats open access publishing in a canonical way as a game between scientists 
and publishers. The behavior of oscillation on closed orbits is a rather rough 
and simplified view of the whole scenario. We do not consider this model to be 
universally valid, but throughout a certain period of time it seems to allow a 
reasonably satisfying description of insights which are intuitively clear. 

This model can be modified. For instance, taking into account that both the 
parameter r reducing the scientist's reputation in case of open access publishing 
and the value of i expressing the reduced impact in case of non-open access 
publishing may change in the context of the game. For example, if the rate of 
open access publishing scientists increases, r will be supposed to decrease etc. 
Playing a little bit and slightly altering the model presented here may result in 
further refinements. 

Moreover, the huge number of evolutionary game-theoretic approaches in the lit- 
erature allows a lot of modifications and further developments of these concepts. 
They vary from biological game theory to concepts that are more applicable in 
economic contexts. These models are good sources and suitable to establish fur- 
ther improvements in a more fundamental manner, see e.g. [3l[6] and [8] as well 
as the references therein. 
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